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The quantum- statistical cluster expansion method of Lee and Yang is extended to investigate off- 
diagonal long-range order (ODLRO) in one- and multi-component mixtures of bosons or fermions. 
Our formulation is applicable to both a uniform system and a trapped system without local-density 
approximation and allows systematic expansions of one- and multi-particle reduced density matrices 
in terms of cluster functions which are defined for the same system with Boltzmann statistics. Each 
term in this expansion can be associated with a Lee- Yang graph. We elucidate a physical meaning of 
each Lee- Yang graph; in particular, for a mixture of ultracold atoms and bound dimers, an infinite 
sum of the ladder-type Lee- Yang 0-graphs is shown to lead to Bose-Einstein condensation of dimers 
below the critical temperature. In the case of Bose statistics, an infinite series of Lee- Yang 1-graphs 
is shown to converge and gives the criteria of ODLRO at the one-particle level. Applications to a 
dilute Bose system of hard spheres are also made. In the case of Fermi statistics, an infinite series of 
Lee- Yang 2-graphs is shown to converge and gives the criteria of ODLRO at the two-particle level. 
Applications to a two-component Fermi gas in the tightly bound limit are also made. 
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I. INTRODUCTION 

In 1957, Lee and Yang showed that the thermodynamic 
functions for quantum-mechanical systems of particles 
obeying Bose-Einstein or Fermi-Dirac statistics can be 
expressed in terms of cluster functions for the same sys- 
tems obeying Boltzmann statistics (Tj-Q • This method is 
based on a cumulant expansion of the grand partition 
function, and hence referred to as the Lee- Yang clus- 
ter expansion method. They subsequently applied this 
method to a dilute system of hard-sphere bosons and dis- 
cussed the A- transit ion of liquid helium-4 [3, 4]. However, 
since a cluster expansion is equivalent to a virial expan- 
sion, it is more suitable for dilute quantum many-body 
systems than dense systems such as liquid helium-4. Un- 
fortunately, back in those days there was no such dilute 
quantum many-body system. Since then, the Lee- Yang 
cluster expansion method has been almost completely 
forgotten, with only a few exceptions (5|-[Tq|. Now, how- 
ever, it seems that the method is best suited to describe 
ultracold dilute atomic gases. In the present paper, we 
extend the Lee- Yang cluster expansion method, which 
was originally developed to describe a system of hard- 
core bosons and fermions with small scattering lengths, 
to describe ultracold atomic gases of bosons and fermions 
with large scattering lengths. 

The Lee- Yang cluster expansion method is a cumu- 
lant expansion of the grand partition function in terms 
of fugacity z := exp(/3/i), and enables us to treat a quan- 
tum many-body problem systematically, where ji is the 
chemical potential and P := l/{k^T) is the inverse tem- 
perature with the Boltzmann constant . Although we 
usually apply the Matsubara Green's function method or 
the Feynman graphs to tackle quantum many-body prob- 



lems, the ranges of validity of these approximations are 
often not well controlled. On the other hand, the Lee- 
Yang method allows us to make a perturbation expansion 
in terms of the number of the interacting atoms, and uti- 
lize the fugacity expansion of the grand potential density, 
which is equivalent to a virial expansion. Moreover, the 
cluster expansion (or virial expansion) for the unitary 
Fermi gas is currently of great interest [liHlJ, because 
the cluster integrals (or virial coefficients) 6^^^ and h^^^ 
have recently been measured using a mixture of ^Li in 
the two lowest hyperfine states [l5[[i6|. Here, the virial 

expansion is ppX^ = 2z ^ ^z'^ + b^^h^ + b^^h^ H , 

with the pressure p and the thermal de Broglie length 
A := (27rfi^/3/m)-^/^. The primary purpose of this paper 
is to identify the criteria of the off-diagonal long-range 
order (ODLRO) of the reduced density matrices in Bose 
and Fermi systems |17l, il8|] based on the Lee- Yang cluster 
expansion method. 

This paper is organized as follows. In Sec. II, the quan- 
tum cluster expansion method is extended for the case 
of a two-component Bose or Fermi system of interact- 
ing particles. We represent the grand partition function 
and the one- and two-particle reduced density matrices 
in terms of cluster functions obeying Bose-Einstein or 
Fermi-Dirac statistics. Our extension is formally appli- 
cable to a uniform system and a trapped system without 
local-density approximation. The derivation and mathe- 
matical justification of our method is given in Appendix 
[Al In Sec. HI, we review the rules by which quantum- 
statistical cluster functions can be computed from the 
knowledge of quantum cluster functions obeying Boltz- 
mann 

classical- (i.e., Boltzmann-) statistical cluster func- 
tions. This formulation was first developed by Lee and 
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Yang for the case of a one- component system Q and 
it was extended by Pathria and Kawatra for a multi- 
component system In Sec. IV, we show that the 
logarithm of the grand partition function as well as the 
one- and two-particle reduced density matrices can be 
expressed as a sum over a particular set of Lee- Yang 
graphs. For the case of a one-component uniform sys- 
tem, the expressions of the grand partition function and 
the one-particle reduced density matrix were obtained by 
Lee and Yang Q , and the two-particle density matrix was 
obtained by de Dominicis [6[. We generalize these results 
to the case of multi-component uniform and trapped sys- 
tems. In Sec. V, we elucidate a physical meaning of Lee- 
Yang primary graphs, in particular, that of an infinite 
series of the ladder- type Lee- Yang primary graphs, the 
sum of which leads to Bose-Einstein condensation (BEC) 
of tightly bound diatomic molecules below a critical tem- 
perature [9]. In Sec. VI, we give the criteria of ODLRO 
based on the Lee- Yang cluster expansion method for the 
case of Bose-Einstein statistics. Applications to a Bose 
gas of hard spheres are also made. In Sec. VII, the cor- 
responding criteria for Fermi-Dirac statistics are given. 
Applications to a two-component Fermi gas in the tightly 
bound limit are also made. In Sec. VIII, we summarize 
the main results of this paper. The detailed proofs of 
several formula are relegated to the appendices to avoid 
digressing from the main subject. The relations between 
the Lee- Yang and our notations are hsted in Appendix 

o 



II. FORMULATION FOR A 
TWO-COMPONENT QUANTUM SYSTEM 

A. Grand partition function 

We consider a two-component gas of bosons or 
fermions whose components are either internal states 
or different atomic species, and distinguished by indices 
cr =1, |. A generalization to a general multi-component 
system is straightforward [see, for example, Ref. Q]- To 
make the wave function square integrable, a system is 
considered to be in a uniform finite volume V = or 
in a trap with an effective volume '^trap- In the case of a 
uniform system, the Hamiltonian is given by 



J2 (K-'^i.'l)^ 



(1) 



where rria- and No- are the atomic mass and the particle 
number in state a. Here we set h = 1. In the case of 
a trapped system, the Hamiltonian should include addi- 
tional terms Xli^ ^trap(rij + Xli^ 

(r^^). We assume 

that the system is kept at constant inverse temperature 
P = (/cb^)~^, and that the system is allowed to exchange 



two kinds of atoms with a reservoir at given fugacities 
and , where the fugacity {cr ="^^1) is related to the 
chemical potential per atom /i^r by z^ = exp(/3/icr)- 
We introduce the function 

:= (iV^!)(7V;!) ^ V^,(l^,...,Ar|;l^,...,7V|) 



x^,m,...,NrAu...,N^)e-^^% (2) 

where a := 1^ := ri^ := {ri^^x,ri^,y,ri^,z), 

and 1^ := vy^ := (ri/^,^, r^^,^, ri/^,^), etc., and ipi 
and Ei are the normalized eigenfunction and the cor- 
responding eigenvalue of H^^^'^^\ respectively. The 
summation in Eq. ([2]) extends over all eigenfunctions 

t/ji G H^^^^ ^ ^a[^\ where the subscript = S or 

A indicates that the A^cr-particle Hilbert space 1-6^/^ is 
symmetrized or antisymmetrized. 

To obtain the logarithm of the grand partition func- 
tion, we define the matrix elements of cluster functions 
^{N^,Ni) coordinate representation as follows [19i] : 

(ll|Wi°-i)|l;)^(l^|l7l°>i)|l;), 
{4,2^|W-(2'°)|l,,2t)^{4,2^|f/^°)|lt,2t) 

+ (4|[/^°)|l,)(ll|C/r)|l,), 
(ll,2^|W^i°''^IU,2;)^{li,2^|f7l°>2)|l,,2;) 

+ {l'M°''^\h) (2l|C/^i)|2;), 

+ {4|C/(i-°)|lt>(2^,3^|f/f°)|2t,3t) 

+ {2^|C/(i'°)|2t>(3^,4|f/i''°)|3t,lt) 
+ (3^|/7(i'°)|3t)(l^,2^|[/(2-°)|1^2t) 
+ {4|C/(i'°)|lt>(2^|C/(i'°)|2t>(3^|C/(i'°)|3t),etc. 

(3) 

The grand partition function is given by 



N^l Nil 



Nf=0 N^=0 

X m,...,Nr,li,...,Ni\ 



(4) 



Wr^'^^^|l^...,iVt;l;,...,iV,), 



where J£^^'r^ := J(fri^---Jd^rNt and Jd^^ri := 
We can show that 



oo oo 



(5) 
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where cluster integrals are 
^(O'O) := 0, 



(6) 



mm 

The proof of Eq. ^ is given in in Appendix [All. 

Once the cluster integrals are obtained, we can evalu- 
ate thermodynamic quantities. For the case of a uniform 
system in a finite volume V = L^, the equilibrium pres- 
sure p and the particle- number density pa- := (fia) of the 
system are given by 



^ ^ CX) CX) 



13 v^oo V 



1-^=0 1^=0 



and 



^ CO CO 



(8) 



Similarly, for the case of a trapped system, the particle 
number No- of the system is given by 



oo oo 



(9) 



B. One- and two-particle reduced density matrices 

In terms of second-quantized field operators, the one- 
and two-particle reduced density matrices can be ex- 
pressed as 



(rcT|pW|rV) := (*t (r) (r')), 



(riCTi , r2(J2 \r'i(T[ , rjCTa) 
:= (*t^(ri)*t.(r2)*.^(ri)*.i(r;)), 



(10) 



(11) 



where {O) is the expectation value of O over the grand 
canonical ensemble defined by 



(12) 



and ^J-(r) and ^cr(r) are the creation and annihilation 
operators. Here, H and N^j := J (i^r^J.(r) ^cr(r) are, re- 
spectively, the second-quantized operators for the Hamil- 
tonian and the total number of particles with spin a. 
[See, for example, Ref. [20|, in particular, Eqs. (2.1.6) 
and (2.4.2)] 



The one-particle reduced density matrix can be shown 
be related to U^^'^^^ by 



X (y,l^,...,/^ - 1;1;,...,/;| 



(13) 



and 



(xibW|yi> 

CO CO 



^ ^ - I "4- ^ ^ 



d3(U-l) 



(14) 



c/^''^^|it,...,;t;x,i;,...,/i-i>. 



The proofs of Eqs. (fT3|) and (fT4|) are given in Appendix 
[32. Clearly, (x t |y 1) = (x ; |y t) = 0. 

The two-particle reduced density matrix can also be 
shown to be related to Ua^'^^^ by 

(xi t,x2tl/o('^|yi t,y2 1) 

= {xitb''^|yi t>(x2tb^'^|y2 1> 



oo oo h h 



EE 



ih - 2)! h 



!U /^'^''"'^rt/rf''^r, (15) 



X (yi,y2,lt>--->^t -2;l4,,...,^4,| 

?7a'^'''''''|xi,X2, 1^, . . . ,1^ - 2;li, . . . , li), 

(xi i,x2 i 1/5^^^ |yi i,y2 i) 

= (xa|p^'^|yi i)(x2i|pW|y2 i) 

oo oo p p 

/7^^^'^^^|lt,...,/t;xi,X2,l;,...,/;-2), 

and 

(xi t,x2;|p(')|yi t,y2 i) 
= e(xi t,x2;|p(')|y2 ;,yi t) 
= e(x2 i,xit|p('Vi t,y2 i) 
= (X2 ;,xit|p('V2 ;,yi t) 
= (xit|p^'^|yi t)(x2;|pW|y2 i) 

oo oo h p p 



EE 



u^i (^t - 1)! - 1) 



X (yi,l^,.. .,1^ - I;y2,l4., . .. ,^4. - 1| 

f/i'^''^Vi,it,---,^T-i;x2,u,...,Zi-i) 



(17) 
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where e = +1 for a = S (Bose-Einstein statistics) and 
6 = —1 for a = A (Fermi-Dirac statistics). The proofs of 
Eqs. (p!5]) - (pT|) are given in Appendix [Al3. Clearly, the 
other components such as (xi t7X2 i Ip^^^Iyi ijy2 i) 
vanish. 

A generalization to an TV-particle reduced density ma- 
trix is straightforward. For example, the three-particle 
reduced density matrix is 

(xi t, X2 1, X3 ; \p^^^ |yi t, y2 1, ys ;) 

= (xi t IP^'^Iyi t)(x2 1 \p^'^\y2 t)(x3 ; \p^'^\ys i) 

+ (xi t Ip^'Vi t)(x2 t,x3 ; |p('V2 t,y3 ;) 

+ (X2 1 |p^'^|y2 t)(xi t,x3 ; |p(^)|yi t,y3 1) 
+ (X3 ; |p(^)|y3 ;)(xi t,x2 1 IP^'^Iyi t,y2 1) 

X (yi,y2,lt,---,^t-2;y3,li,---,^;-l| (18) 

V^''''^ |xi, X2, It, . . . , /t - 2; X3, 1;, . . . , - 1). 

The formulation given in this section is expressed in 
the form of an infinite series. Thus, one might question 
the mathematical rigor of the formulation. Actually, the 
derivations of the formulation can be mathematically jus- 
tified, as shown in Appendix [Al-4. 

III. U^^'^^^ IN TERMS OF U^^^'h^ 

In the preceding section, we derived the formulas of the 
grand partition function and the one- and two-particle re- 
duced density matrices. To evaluate these quantities, we 
must calculate Ua^'^'^\ In this section we review the rules 
by which Ua^'^^^ can be computed from the knowledge 
of certain /7^^^'^^^ functions for the same system obeying 
Boltzmann statistics. These rules were first established 
by Lee and Yang for the case of a one-component sys- 
tem jll and was extended by Pathria and Kawatra for a 
multi-component system [7j. 

A. Boltzmann statistics 

For Boltzmann statistics, we introduce the function 

(1^, . . . , iv;; 1^, . . . , iv||M/(^t.^^) lit, . . . , iVt; • • • , ^^i) 

xi;im,...,Nr,li,---,Nl)e-^''\ (19) 

The summation on the right-hand side of Eq. (p!9|) runs 
over all eigenfunctions ipi G T-L^^^^ ®1-L^^^\ where l-L^^""^ 
is the A^o-particle Hilbert space. 



Similarly to Eq. (j3j), we define the matrix elements of 
cluster functions U^^^'^^^ in the coordinate representa- 
tion as 

(i^iw(i-°)|it) = (4|[/(i-°)|it), 

(l^|W-(°'i)|lt) = {ll|f/(°-i)|lt), 
(4,2^|W(2'0)|lt,2t) = (4,2^|[/(2'°)|lt,2t) 

+ {4|[/(i-°)|lt) (2^|C/(^'0)|2t), 
{4;ll|W-(i-i)|lt;lt) = {4;ll|[/(^-^)|lt;U) 

+ {4|[/(i-°)|lt)(ll|C/(°-i)|lt), 
(l^,2^|M/(0-2)|lt,2t) = (l^,2^|[/(°-2)|lt,2t) 

+ (l^|[/(°-i)|lt) (2^|C/(°'i)|2t), 
(l^,2^,3^|Vr(3.o)|lt,2t,3t) 
= (l^,2^,3^|[/(3-0)|l^,2t,3t) 

+(l^|/7(i-°)|lt)(2^,3^|/7(2.o)|2t,3t) 
+(2^|/7(i'°)|2t)(3^,4|C/(^'°)|3t,lt) 
+(3^|C/(i'°)|3t){4,2^|C/(2'°)|lt,2t) 
+(4|/7(i'°)|lt)(2;.|C/(i'°)|2t)(3^|[/(i'°)|3t),etc. 

(20) 

B. f/^'^''^* in terms of U'^h'^i) 

We now formulate the rules by which Ua^''''^^ can be 
computed from the knowledge of U^^^''-^\ Such rules re- 
sult from the fact that Ua^'^'''^ and {/''f''^^ are defined 
in terms of wl'^''"'''' and VT^t.^i)^ respectively, which, in 
turn, are related to each other through 

...,N[;ll,..., iVjlT^r^'^^Vt, .■.,Nr,U,..., TVt) 

= E E ^''^'^(^(it). • • • ' ^(^;); Q(il)' • • • ' Q(^I)I 

W(^t,A'.)|i^,...,iV^;l^,...,iV^), (21) 

where e = +1 for a = S (Bose-Einstein statistics) and 
e = — 1 for a = A (Fermi-Dirac statistics). Here, P and 
Q denote permutations, and (—1)^ and (—1)^ take on 
1 or — 1 for even or odd permutations. Equation ([2T]) is 
proved in Appendix [Bl It can be seen from the structure 
of Eqs. ([20]) and ^ that ui^^'^^^ may, in general, be 
expressed in terms of the sum over connected products 
of Ta^'^^^ functions, where 

/7(^t,U)|i^,...,/^;l^,...,/^). (22) 
Here, we only show the results: 



5 



Rule — To calculate (r'^^, . . . , rj^; r'^^, . . . , r[^\Ua^'^'^^ 
|ri^, . . . , Ti^; ri;, . . . , r^), we first consider a grouping of 
the + li integers 1^, • • • , 1;, • • • , h- 

{K)(^t) • • • }{(ct^^t)(et/t) • • • }{(^t^tH) •..}..• 

{Uth)ih^i) ' ' • }{K^m) • • • }{(^tnn) ...}... 

(23) 

where a^, 6^, . . . , g^, . . . and k^, m;, . . . , . . . are per- 
mutations of these and integers, respectively. In each 
pair of the curly brackets, there are mn^,n^(= 0, 1, 2, ... ) 
round brackets with up-spin integers and down- 
spin integers which are subject to 



n^=0 n^=0 



(24) 



Within each pair of the round brackets in (|23|) . the inte- 
gers for the same spin are arranged in ascending order, 
i.e., c < e < g < h < o < etc., and within 
each pair of the curly brackets, the round brackets are 
arranged in such a manner that their first entries are 
placed in ascending order. We then form the sum 



E({(rAtin^'°^|rat>(r'Btin^'°^M---} 

|T(2-°)|r,t,r.t)--- ••• 

r . (25) 

x{(r'^,;r'^JTii-i)|r,t;rfe,>...}... 

x{(rH|Tr^|ru){rW|Tr^^|r«4>---}-' 

where A, 5, . . . , J, . . . and i^, . . . , T, /7, . . . are permuta- 
tions of 1, 2, . . . , and 1, 2, . . . , Z^, respectively. Because 

we use Ta^'^^^ in Eq. ([22]) . only one permutation will be 



included in the sum in Eq. (|25]) among all the permuta- 
tions A, 5, . . . , J, . . . which differ from each other only 
in the relative positions of numbers within the same bra 
(e.g., (rj^.^,r^.^| and (r^^,rj^^| ). The sum in Eq. ([25]) 
then runs over all permutations, provided that upon set- 
ting = r^^ and r-^ = r^^ for all z, the summand in 
Eq. ([25]) cannot be written as a product of two factors, 
one of which depends only on some, but not all, of the 
coordinates ri^, . . . , r^^ and ri^, . . . , r^, while the other 
depends only on the rest of these coordinates. The sum 
of all expressions in Eq. ([25]) over the different groupings 



in Eq. 



then gives Ua 



We give some examples. 



Example 1: 



Example 2: 

(l^,2^|[/(2.o)|i^^2t) = (l^,2^|Ti2.o)|i,,2^) 

+ e(2^ini-°)|l,>(4ini-°)|2,). 

Example 3: 

(4;ll|/7(i-i)|lt;l;) = (4;ll|Tii.i)|lt;l;). (28) 
Example 4 ' 

(4,2^;3i|l7^i)|lt,2t;3;) = {l\,2\;S[\T^^''^^,2r,h) 
+ (2^|Tii.°)|lt)(4,3^|TLi-i)|2t,3,) 

+ {4|T(i'°)|2,){2i,3^ini-i)|l,,3,). 

(29) 

For reasons which will become clear later, it is useful 
to define a function T^^ related to T^^'^^^ as follows. 
First, we define (f i, . . . , fz) := (ri^, . . . , r^^, ri;, . . . , r/;), 
(f'i,...,fj) := (r^^^, . . . , rj^, r'^^, . . . , rj^), ai = ••• = 
ai^ =t and a^^+i = • • • = a/ =|, where I := l^ ^ 

Then, we define a function T^^ as 

(^Q(i)'^Q(i)5---'^'q(0'^Q(0I 

:=6^6^(f;,...,f;|TL^^'^^^|fi,. ..,£•,) (30) 
= e^e^{r[^, . . . , rj^; r'^^, . . . , 

xi^^'^^ Vit' • • • ' '^^t; ^u, . . . , r^), 

where P and Q are any one of /! permutations. 
We give a few examples. 



Example 1 : 



(r',t|TW|r,t) = {r'|TLi-°)|r>. (31) 
(r',;|TW|r,i) = {r'|n°-i)|r). (32) 



Example 2: 



Example 3: 



(l^|C/^°)|l^) = (4|Tii.°)|lt). 



(26) 



{r;,t;ri,t|Ti2Vi,t;r2,t) 

= e(r^,t;r'i,t|TL2)|ri,t;r2,t) 

= e(r;,t;r^,t|TL2)|r2,t;ri,t) (33) 

= (r^,t;r'i,t|Ti2)|r2,t;ri,t) 

= (r;,ri|Ti2.0)|ri,r2). 

{r'i,t;r^,;|Ti2)|ri,t;r2,;) 

= e(r^4;r'i,t|TL2)|ri,t;r24) 

= e(r;,t;r^4|TL2)|r24;ri,t) (34) 

= (ri4;r'i,t|TL'^|r2,;;ri,t) 

= (r;;r^|Tii-i)|ri;r2). 
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IV. GRAPHICAL REPRESENTATION 
A. Primary ^-graphs 

Using Eqs. (|5]), (|6]), (US])- ([III) and the rules introduced 
in the preceding section, the grand partition function and 
the one- and two-particle reduced density matrices can 
be expressed in terms of the sum over connected products 
of Tq, functions. The exact character of this sum is most 
simply described in terms of primary ("-graphs introduced 
by Lee and Yang Q for the grand partition function and 
the one-particle reduced density matrix, and by de Do- 
minisis [6| for the two-particle reduced density matrix. A 
primary ("-graph is defined as follows: 

Definition. — A primary ^-graph {( = 0, 1, 2, ... ) is a 
graphical structure which consists of a collection of ver- 
tices connected by directed lines, with (" external incom- 
ing lines and ( external outgoing lines. Here, a line that 
has vertices at both ends is called an internal line; oth- 
erwise, it is called an external line. All external lines 
are considered distinguishable. Each vertex, called the 
/-vertex (/ = 2,3,...), connects / incoming lines and / 
outgoing lines. A primary (^-graph must include at least 
one vertex and one line, and all parts must be connected 
(i.e., there must be a path from any one vertex to any 
other vertex). Two primary graphs are different if their 
topological structures are different. 

To each of these graphs we assign a term which is de- 
termined by the following procedures: 

(i) Associate with each internal line a different integer 
i {i = 1, . . . , A^) and the corresponding coordinate 
and spin (r^,<j^). Associate with each external line 
some prescribed coordinate and spin. 

(ii) Assign to each /-vertex, a factor 
Bi B2 . . . Bi 

=4^4^(5i,...,5,|T(^)|Ai,...,A,), (35) 
A1A2...A1 

where Ai and Bi represent the coordinates and 
spins associated with its incoming and outgoing ith 
lines (i = 1, . . . , /), respectively. The number of up- 
spins associated with its incoming lines is the same 
as that of outgoing lines, and it is denoted as l^. 
The same is true for down-spins Z^, with / := 

(iii) Assign a factor 1/S to the entire graph, where S is 
the symmetry number and is defined as follows: 

Consider all N\ permutations of the positions of N 
integers associated with the lines. The total num- 
ber of permutations that leave the graph topologi- 
cally unchanged gives the symmetry number of the 
graph. The symmetry numbers are listed under 
each graph in Figs. [TJ [21 etc. 



(v) Assign a factor —1 to the entire graph, if the per- 
mutation 

Ai ^ Bu A2^B2, A^^B^, ... (36) 

from all the initial coordinates into all the final co- 
ordinates of all the vertex function T^^ taken to- 
gether is odd. 

The term that corresponds to each graph is given by 

at^,.. J J (37) 

X [ product of all factors in (ii) — (v) ] . 

B. Primary 0- graphs 

In terms of these primary graphs, we can write the 
grand partition function [3] as 

logS (/3, z) = [all different primary 0-graphs] , (38) 

to which each graph contributes a term given by Eq. ([37]) . 
Equation ([38]) is illustrated in Fig. [H We can write 
Eq. ([38|) explicitly as 

logS(/3,^) 



= z^Tr ([/(I'O)) + -i Tr (U^^^'^^f + Tr (U^^^'^^f + • • • 
2 3 
2 3 




+ ••• , 

where TrX := /d^r (r|X|r), TrX^ := Jd^n J(fr2 
(ri|X|r2){r2|X|ri>, etc. 



C. Primary 1- and 2-graphs 

It is useful to define functions {r' a' \r]^^^ \ra) and 
(r'^cr^, r2cr2|7^^^^ |ri<Ji, r2cr2) related to the one- and two- 
particle reduced density matrices by 

(rV'|7/(i) |ra) := 5{r' - r)5^^. + e{ra\p^'^ |rVO, (40) 

and 

{l'2'|r?(2)|l2) 

-(l'|r?«|l>(2'|r?«|2)-e(2'|r?«|l)(l'|r?«|2) 
= (l'2'|p(2)|i2) (41) 

- (l'|pW|l)(2'|pW|2) - e(2'|pW|l)(l'|pW|2), 
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(a) 
log 3 



(1) (2) (3) 



+ 





+ ••• 



+ 



(2) 



= tr(log /7o) + X )+ Ok 



(b) 



(2) 

+ i 4 I ; + ••• 

(2) 



(6M 



o 

+ i + o + 
o 



(1) (1) 



(1) 



FIG. 1: (a) Expression of logS as the sum over all differ- 
ent primary O-graphs (solid curves) and as the sum over all 
different contracted O-graphs (dotted curves). The symme- 
try number is shown under each graph, (b) Expression of 
(rV^|?7Q"'^'* |rcr) as the sum over different primary 1-graphs. 



where 1,2 and l',2' represent both the coordinate and 
spin of the particle such as 1 := (ri, <Ji) and 1' := (r^,^-^). 

In terms of second-quantized field operators, rj'^^^ and 
T]^'^^ can be written as 



(rV|7?«|m) = (*,(r')*t(r)), 



(ri<T,r2<T'|r?(2)|r;a,r^<7') 



(42) 



(43) 



These functions can be shown to be related to 1- and 
2-graphs 0, by 



{r'a'\rj(^)\ra) = 6{r - r)(5^^/+e ^ [all different 

primary 1-graphs] , 



(44) 



and 

(l'2'|r?(2)|l2> 

- {l'|r?W|l){2'|r?W|2) - e{2'|r?(i)|l)(l'|r?(i)|2> (45) 
= [all different primary 2-graphs] . 

D. Contracted (^-graphs 

In Fig. [1] (a), we have arranged the sum so that the 
different primary graphs in the same row have the same 
number of /-vertices for / > 2 and they differ only in the 
total number and relative positions of 1-vertices. We see 
that all the terms in the same row in Eq. ([39]) can easily 
be summed up. Thus, Eq. (j39j) becomes 

logS = Tr [log (7?^'^) 
2 

x(rithr|ri tXrs theirs t) 
x{rit|r?riri t>(r2 i |r?rir2 i) 

+ ■■■ , 

where 

{r't|r?^'^|rt) := {r'|(/ - e^C/(i'°))-» 

= 5(r'-r) + ez^(r'|C/(i.o)|r) (47) 

+ 4y"dV{r'|C/(i'°)|ri)(ri|C/(i-°)|r> + ..., 

(r'i|r?W|r;) := {r'\{I - ezU^°''Y'\r) 

= 5(r'-r) + ez;(r'|[/(°-i)|r> (43) 

+ zl I d3ri(r'|C7(°'i)|ri)(ri|/7(°'i)|r) + • • • . 

It is therefore convenient to represent the primary graphs 
in the same row as in Fig. [1] (a) by a single structure 
which is obtained from any one of these primary graphs 
by simply eliminating all 1-vertices (see Fig. 1 (b)). 

We now introduce a contracted graph. A contracted 
graph (or, contracted ("-graph) has the same topological 
structure as a primary graph (or, primary ("-graph) ex- 
cept that a contracted graph does not have any 1-vertex. 
To each contracted graph, we assign a term which is de- 
termined by the same procedures (i) -(v) and the follow- 
ing additional rule: 

(vi) Assign a factor t^q^^ to the zth internal line. 
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In terms of these graphs, Eq. (|46|) reduces to 



(a) 



logS(/?,^)=Tr[log (r/W) 



[ah different contracted 0-graphs] 



(49) 



This sum is also ihustrated in Fig. [T] (b) in which sohd 
hues represent the primary graphs and dotted hues rep- 
resent the contracted graphs. 

Similarly, Eqs. (j44]) and (j45j) become 

+ e [all different contracted 1-graphs] , (50) 

and 

(l'2'|r?(2)|l2) 

- (l'|r?(i)|l){2'|r?W|2) - 6{2'|r?W|l){l'|r?W|2) (51) 
= [all different contracted 2-graphs] . 

Equation (|5Q|) is illustrated in Fig. [2l 

For further discussions as well as examples of primary 
and contracted ^-graphs, the reader is referred to Sec. 3 
of Ref. ^] and Sec. II - (ii) of Ref. 



E. One-component system 

In this subsection, we consider a one-component Bose 
or Fermi gas. We can establish the formulae for the one- 
component gas by substituting 2:^ = in the correspond- 
ing formulae of the two-component gas. From Eqs. ([39|) 
and (|36|), the grand partition function is 

logS(/3,z) 

= z Tr ([/(I)) + -z^ Tr [V^^^f + -z'' Tr (IJ^^^f + • • • 
2 3 



Jdhidh2 (ri,r2|Ti2)|ri,r2) 
X ('i^2^z3(ri|/7(i)|ri) + --. 



= Tr[log {rii'^y 

+ i Jdh.d'r, (ri,r2|n2)|ri,r2) 

X (ri|7?^^Vi){r2|??o^V2) 
+ ■■• , 

where z := z-f C/C) := tL'^ := tL''°^ and 

(r'|,?W|r) :=(r'|(J-e^;7(i>°))-». 



(52) 



i 


A ....... / 
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A ■' / 
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^ ■■• 
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•••■/k '-"^-' ■ 
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(1) 
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+ !>. 


+ A.-K 


-.. + •• 






y 










■ (1) 




+ ■■■ 







y + ••• 



y + 



(b) 




FIG. 2: (a) Expression of (rV'lr?*^' |m> as the sum over all 
different contracted l-graphs, where the sum is rewritten in 
the form of a geometric series (last line), (b) Expression of 
(rV'l^^^^cr) as the sum over all different irreducible con- 
tracted 1-graphs. The corresponding symmetry numbers are 
shown under the graphs. 



V. PHYSICAL MEANING OF PRIMARY 
GRAPHS 

For a uniform system in a finite volume V, the cluster 
functions in momentum representation are 



(k;,...,k;ic7('t^)|ki,...,kz) 



1 ' 



(53) 



X {r;,...,r;|C/('t.U)|ri,...,r,), 
where (ki,...,ki) = (kj, . . . , k;.^; k;.^+i, . . . , k;), etc. 



(2) (2) 



1^2 



1^2 



interaction quantum exchange 

between 2 particles between 2 particles 

FIG. 3: The second-order terms in fugacity z for the primary 
0-graphs [see Eq. (|6T]) ]. 



For later use, we evaluate the cluster functions U^^^ for 
a one- component system and /7^^'^^ for a two-component 
system. We can evaluate IJ^^^ for the one-component 
system as 



(55) 



Next, we evaluate [/^^'^^ for the two-component system. 
The two-particle Hamiltonian is 

if (1.1) ^ ^k2+— ki+« = - — i — ,k2+f(-'), 



2mt 2m i 



^(rel) 



2m 



2(m^ + m;) 



(56) 
(57) 



where we introduce the center-of-mass momentum K := 
ki + k2, the relative momentum k := (ki — k2)/2, and 
reduced mass m* := m^m^/ (m^ + m^). The cluster func- 
tion is 



(k;k'2|[/(i>i)|kik2) 

Stt^ 



V 



<5k,K' e-''^ /{2(r™,+r™,)} (k'|[/(rel)|k)^ 



(58) 



where 



(k'|/7(-')|k) := (k^-z'^'^'lk)- Vk'e-'"^'/^^™-)^ (59) 

The function (k'|/7^^®^^ |k), which is independent of the 
volume, includes the effect of the interaction and can 
be calculated by the eigenfunction and the eigenvalue of 

^(rel)^ 

In the following, we list four simple examples for a 
uniform and equal-mass (m^ = m^ = m) system in a 
finite volume V. 

Example 1 (Quantum exchange and interaction): 

The grand partition function of a one-component Bose 
or Fermi gas is expressed as a sum of primary 0-graphs. 
Let us consider any one of the primary 0-graphs. It 
has several /-vertices and their orders are denoted as 
/i, . . . , Z^. The total number of the orders, X^ILi 
responds to the order of fugacity z of the corresponding 



n=l 

iri) 



(1) (2) (3) 

k 



+ 



(2) 
ko ki 



ko ki ko 



k ^ ^ ^2 ^ ^ ^ki^, kg ^ ^ ^ 

k ki ko ki ko ko 



+ 



quantum exchange quantum exchange 
between 2 particles between 3 particles 

FIG. 4: Sum over different primary 0-graphs which include 
the 1- vertex only [see Eq. (|63|) ]. 



primary 0-graph. In general a contribution from N par- 
ticles is expressed as the A/'th-order term in z. Hence, the 
primary 0-graph represents the many-body effect among 
Xir=i particles. The one-particle cluster function is 
given by 



-^A-3exp[-(r'-r)V(4/3)], (60) 



{r'|C/(i)| 



where A := (27r/3/m)^/^ is the thermal de Broglie length. 
This expression can be better understood if we separate 
the effect of the interaction from the quantum-exchange 
effect. To provide a concrete example, we consider the 
second-order in fugacity z for the primary 0-graph [see 
Fig. [3]: 

^ Jd'r.d'r, ((ri,r2|T(^Vi,r2) 

+ e(r2|[/^^Vi)(ri|t/^'V2)). (61) 
Here, the first term in the integrand is 



(ri,r2|Ti2)|ri,r2> 

= (n, r2|C/(2) ^2) + e(r2, ri |f/(2) ^^^^^^^^ 



(62) 



where the cluster functions /7^^^ contain only the 
effect of the interaction, because they are defined with 
unsymmetrized wave functions. The second term in the 
integrand of (|6T]) describes the particle-exchange effect. 
According to Eq. (|6Q|) , the relevant scale of the term 
e(r2|/7(^)|ri)(ri|[/(^)|r2) is the thermal length A. We 
can thus divide the effect of the two-particle correlation 
into the quantum-exchange effect of the Bose-Einstein 
or Fermi-Dirac statistics and the effect of the interaction 
with Boltzmann statistics. 

Example 2 (Ideal Bose gas and BEC): 

For an ideal Bose gas, BEC is caused by the quantum- 
exchange effect such as the first term in Fig.[3l For nonin- 
teracting particles, it is clear from Eq. (|2Q|) that U^'^^ = 0, 
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U^^^ = 0, etc. Therefore, we obtain from Eq. (j38]) 
log E = ^ [all different primary O-graphs] 
= ^^(k|C/(i)|k) + f E (k2|f/(^)|ki)(ki|C/(i)|k2) 



ki,k2 



J2 {k2|f/(^)|ki>(k3|C/(l)|k2>(ki|C/(l)|k3) 



ki,k2,k3 



(63) 



where each term in the sum corresponds to the primary 0- 
graph in the sum in Fig. HJ The summation extends over 
ki e 27rZ^/L, where we use the symbol Z = {0, ±1, . . . } 
for a set of integers. As mentioned earlier, the order 
(power) of fugacity z indicates the number of the particles 
involved. For example, the first term on the right-hand 
side of Eq. (|63|) corresponds to the free motion of one 
particle, the second term corresponds to the quantum- 
exchange effect between two bosons, and the third term 
corresponds to the quantum- exchange effect among three 
bosons. Substituting Eq. (|55]) in Eq. (|63|) . we obtain 



-n/3k^/(2m) 



k n=l 



(64) 



which agrees with the well-known result for ideal 
bosons. At z = e^^ = 1, the summation in Eq. ([64|) 
has a singularity at k = 0, which is characteristic of a 
Bose-Einstein phase transition. 

Example 3 (Ideal Fermi gas and Fermi surface): 
For ideal fermions, we obtain 

log S = ^ [all different primary O-graphs] 

= zY,{W''^\^) - f E {k2|f/(^)|ki)(ki|C/(i)|k2) 

k ki,k2 

+ ^ E {k2|f/(')|ki)(k3|C/W|k2)(ki|C/W|k3) 



ki,k2,k3 



(65) 



EE 



/3k2/(2m) 



n 



= Ei«s( 



k n=l 

l + ^e-/3kV(2-) 



which agrees with the known result for ideal fermions. 

Example 4 (Two-component Fermi gas in the tightly 
bound limit): 

We assume that the two-body interaction v is finite- 
ranged and supports one bound state corresponding to 



a diatomic molecular state. Let V^b(k) be the bound- 
state solution for the center-of-mass system of two par- 
ticles with binding energy E\) such that i^*^^®^^'?/^b(k) = 
£^b'0b(k). We assume that the binding energy is much 
larger than the thermal energy, i.e., |£^b| ^ te^- For 
the Bose-Einstein distribution of diatomic molecules to 
develop a singularity, we should have fi < — \E\)\ /2. Us- 
ing Eq. ([59]) . we obtain 



(k'|C/(-')|k>~^b(k')^:(k) e'^l 



(66) 



We can show that a general formula for a central poten- 
tial v{r) is given by 



{k;,k'2|Ti^'^)|ki,k2; 

^ g-/3KV(4m)^j^_j^,^^^(k')t/':(k) e'^l^^l. 



(67) 



The factor e~^^^/^^^^(5K,K' is associated with a free mo- 
tion of the center of mass of molecules. Then, if we con- 
sider an infinite series of primary O-graphs that corre- 
spond to the quantum exchange of molecules, the sum of 
the infinite series leads to the BEC of diatomic molecules 
below the critical temperature. We consider the sum of 
ladder-type primary O-graphs as shown in Fig. [5l The 
number in the parentheses under each graph in Fig. [5] is 
the symmetry number of the corresponding ladder-type 
primary 0-graph. The explicit algebraic expression of the 
sum of the graphs is given as 



n=l ki,...,k2n Cri,...,Cr2ri 



(68) 



(kiCri,k2Cr2|T^^^|k3Cr3,k4Cr4) • • • 
X (k2n- 1 Cr2n- 1 , k2n Cr2n \ T ^ ^ | ki ai , k2 (72 ) . 

Substituting Eq. ([67]) into Eq. (j68]), we have 

^'ad = EE ^(-'^"^^'^"''''/^'"'O" (EiV'b(k)i')" 



K n=l 



^log(l-^VI^-le-^^'/(^-)). 

K 



(69) 



Atz = ^ ^ _ /2), the sum in Eq. 1^ 

develops a singularity at K = 0, which is a manifesta- 
tion of the Bose-Einstein phase transition of diatomic 
molecules. 

Note that the sum of ladder-type graphs shown in 
Fig. [5] gives rise to the quantum-exchange effect of 
the Bose-Einstein statistics of diatomic molecules. In 
the Matsubara Green's function method, a bound state 
emerges as a result of an infinite summation of ladder- 
type Feynmann diagrams which is interpreted as describ- 
ing an effect of repeated binary collisions. In the Lee- 
Yang method, however, since the cluster function t/*^^^ 
is evaluated from a non-perturbative bound-state solu- 
tion, the single cluster function [/^^^ includes the effect 
of repeated binary collisions. In this method, we should 
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oo 
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1 2 




^ 1,2 / \ ^ 1,2,3,4 / \ / \ 
12 12 3 4 

pair pair ^ pair 

quantum exchange between 2 pairs 

3 4 5 6 1 2 



+ 



1 Y< ^y^^ Ny/ 

24,^, A A A 

1 2 3 4 5 6 
pair pair pair 
quantum exchange between 3 pairs 



FIG. 5: Viad [Eq. (|68)) ] as the sum over different ladder- type 
primary 0- graphs. The symmetry number is shown under 
each graph. 



carry out an infinite summation of ladder- type Lee- Yang 
graphs to include a quantum-exchange effect of the Bose- 
Einstein statistics, since the cluster function [/^^^ is de- 
fined in terms of the Boltzmann statistics. 



rest being of the order of 1. If there exists more than 
one eigenvalue of the order of the BEG is said to be 
fragmented 

To establish the criterion of BEG in terms of the Lee- 
Yang cluster expansion method, we first consider the re- 
lationship between the eigenvalues of p*^^^ and r]''^\ It 
is convenient to use the basis set {xi(i*cr)}. The matrix 



elements and r]lj^ are defined by 



(72) 



and 



J2 Jdh|dh'{r'a'\v^'^\ra)x^{ra)x*{r'a'). (73) 



Using Eq. (|iQ]) . we obtain 



(74) 



Therefore, if the system shows BEG, the eigenvalue of 
7^^^^ is of the order of the total number of particles N := 
(TV), i.e.. 



(0 < X, < 1). (75) 



VI. CRITERIA OF ODLRO IN ROSE SYSTEMS 

A. OfF-diagonal long-range order in Bose systems 

The defining characteristic of Bose-Einstein condensa- 
tion is the occurrence of macroscopic occupation of one- 
particle states [13 • The one-particle reduced density ma- 
trix (r(j|p*^-^^ |rV) := (^J.(r)^^/ (r^)) may be expanded in 
terms of its eigenfunctions Xii'^^) with eigenvalues Ui as 

{r<T|p(i)|rV) = ^n.x:(r<T)x.(rV). (70) 

i 

According to Penrose and Onsager [l^, the system ex- 
hibits BEG if one or more of the eigenvalues are of the 
order of the total number of particles N := (TV), i.e., 

m = XiN + o(TV), {0<Xi< 1). (71) 

It can be shown that in an infinite system this condition is 
equivalent to the emergence of an off-diagonal long-range 
order (ODLRO) [TsJ. A single BEG is said to emerge if 
one and only one eigenvalue is of the order of TV, all the 



B. Criterion of BEC in terms of irreducible 
contracted graphs 

We introduce the concept of reducihility of a primary 
1-graph or a contracted 1-graph. 

Definition. — A primary 1-graph or a contracted 
1-graph is reducible if by cutting one of its internal lines, 
the entire graph can be separated into two disconnected 
primary 1-graphs. 

Definition. — An irreducible primary 1-graph or an 
irreducible contracted 1-graph is a primary 1-graph or a 
contracted 1-graph which is not reducible. 

Note that we use the term "irreducible" in a conven- 
tional sense. The same term was used for a different 
meaning by Lee and Yang Q. The correspondence be- 
tween the Lee- Yang's notation and ours is listed in Ap- 
pendix o 

In terms of irreducible contracted 1-graphs, 
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{r^ cr^ \r]^^'> \rcr) can be expressed as a geometric series: 
(rV|7?(i)|m) = (rV|7?^^)|m) 

CTl ,(72 

X {v2cr2\rlo^\va) 
+ Jd^rvd^v,{v'a'\vi'^\ria,){v,ai\K('^V2a2) 

(71 ,...,a4 

X (r2CT2 1??^^^ kacrs) (r3(73 I J^(^^ |r4Cr4) (r4(74 1?7^^^ \ra) 

+ •••, (76) 

where 

= [all different irreducible contracted l-graphs] 
= ^(rV,ricri|T^^^|rcr, ricri)(ricri|7^^^Vi^i) 

ricri 

+ 2;^ ^ (r a ,ricri,r2cr2|T2^Vcr,ricri,r2cr2) 



X (riai|7?^'Vi^i)(r2a2|/?nr2a2) (77) 



This geometric series is illustrated in Fig. O The matrix 
elements rj^^^-- and K^^"^ are defined by 



% :=E/^'r|d3r'(rV'hW|r'^)x.Mx*(r'a'), 



(78) 



and 



^= E j<i'rjd^v'{v'a'\k^%a)Uva)x*{r'a'). 



We can rewrite Eq. ([76|) as 



Then, we obtain 



(79) 
(80) 

(81) 



This equation is illustrated in Fig. [6l Thus, the following 
theorem is established: 

Theorem 1. — 
The following two criteria are equivalent: 

1) The largest eigenvalue of the one-particle reduced 
density matrix p^^^ is 0{N) as ^ 00. 

2) The smallest eigenvalue of the positive definite 

matrix (^r]^^^^ - K^^^ is 0{1/N) as TV ^ 00. 

The former implies ODLRO and the latter can be ob- 
tained from the sum over the contracted graphs. 




FIG. 6: Expression of (rV^|?7*^"'^'* Ira) in terms of an iterated 
equation. 



C. Criterion of BEC in terms of irreducible 
primary graphs 

We can regroup primary l-graphs and estabhsh a cri- 
terion of BEC in a more transparent manner. In terms 
of irreducible primary l-graphs, (rV|7^^^^ |r(j) can be ex- 
pressed as a geometric series: 

ricri 

+ E {rV|/f(i)|ri(Ti)(ri<Ji|i^(iV2<T2){r2<T2|/f(i)|m) 

rir2Cricr2 

+ • • • , (82) 
where 

= ^ [all different irreducible primary l-graphs] (83) 



This geometric series is illustrated in Fig. [71 

Since matrices {r' a' \r]'''^^ \ra) and {r'a'\K'^'^'>\ra) com- 
mute, they can be diagonalized simultaneously. There- 
fore, we can write 

{ra\K(^^\r'a') = E A, x:(ra)x.(rV). (84) 



We can write Eq. ([82]) as 

oo 

n, + l = E(A»)' = (l-Ai)"'. (85) 



Hence, we obtain 



^=0 



A. = 1 



1 



Thus the following theorem is established: 



(86) 



Theorem 2. — 
The following two criteria are equivalent: 
1) The largest eigenvalue of the one-particle reduced 

density matrix p^^^ is 0{N) as A/" ^ oo. 
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FIG. 7: (a) Expression of (rV^|?7*^"'^'* Ira) as the sum over all 
different primary 1-graphs, where the sum is rewritten in the 
form of a geometric series (last line). The symmetry number 
is shown under each graph, (b) Expression of (^v' g'\K^^^\vg) 
as the sum over all different irreducible primary 1-graphs. 



2) The largest eigenvalue of the matrix K^^^ is 
1 + OiXjN) as TV ^ oo. 

The former implies ODLRO and the latter can be ob- 
tained as the sum over the contracted graphs. Thus we 
can use Theorems 1 and 2 to find out whether the Bose 
system exhibits ODLRO based on the Lee- Yang cluster 
expansion method. 



D. Application to a uniform one-component Bose 
gas of hard spheres 

As an example of the formulation given in this section, 
we consider a uniform one-component dilute Bose sys- 
tem of hard spheres of radius a in a gaseous phase and 
evaluate the critical value Zc for the BEG transition in a 
power series of a. Now 

(k'|i^^i)|k) =^(k'|[/(i)|k) 

+ ^(k^ q|T^') |k, q)7^o(q) + 0{c?). (S^) 



Using an explicit form of [/^^\ which was derived in Ref. 
[il, we have 

(k',q|T(')|k,q) = 2(k',ki|C/(2)|j^^j^^^ ^o(^2) 

= _5^^,f^e-''('^?+'^')/(2™) +0(a2). ^^^^ 

From Eq. (jHi) and r?o(q) = (1 - ^e-'^'^'/^^m))-!^ j^^^^g 

^(k,q|T(2)|k,q)r,o(q) = -^Li3(z) I + 0(a2), (89) 
q 

where Li„(z) := -^V^" the polylogarithm. By 

substituting Eqs. (|55|l and (|89)) into Eq. ()87)) . we obtain 



(k'|i4:(i)|k) 



^^^,g-;3kV(2m)^ (l - 4 Li3 {z) ^) + 0(a2). 



(90) 



According to Theorem 2, the critical fugacity Zc satisfies 
the following equation: 

1 = (l - 4Li3 (zc) ^) + 0{c?). (91) 

Hence, we obtain 

'3\ a 



^, = l + 4C( ^ 1 ^+0(a2), 



(92) 



where C(3/2) := Li3/2(1) — 2.612. This result agrees 
with that obtained by Lee and Yang 



VII. CRITERIA OF ODLRO IN FERMI 
SYSTEMS 

A. OfF-diagonal long-range order in Fermi systems 

We consider a system of two-component Fermi parti- 
cles with fixed temperature T and fugacity z. A gen- 
eral criterion of the quantum phase transition of a Fermi 
gas [e.g., superconductivity] is characterized by the oc- 
currence of macroscopic occupation of two-particle states 
[Tsl . The two-particle reduced density matrix 



:= (*|(ri)*|(r2)*4(r'2)*tK)) 



(93) 



may be expanded in terms of its eigenfunctions Xi (1*11*2) 
with eigenvalues n\ as 



(2) 
^ rir2;r'^i 



rx*(ri,r2)x.(r;,ri). 



(94) 



The system shows ODLRO if one or more of the eigenval- 
ues n\ are of the order of the total number of particles 
AT := (iVt) + (iV^), i.e., 



(0 < < 1). 



(95) 
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To establish the criterion of ODLRO in terms of the 
Lee- Yang cluster expansion, we consider the relationship 
between the eigenvalues of p^'^^ and ri^'^\ We define ma- 
trices ri^'^\ and j]^'^®^^ as follows: 



(96) 



r; r2;rir2 



(r;t|p('^|rit)(rUlp('^|r2i), (97) 



and 



The identity matrix I is 

-fr'.r' ;rir2 — -^(r'l " Ti) (5(r2 - Fa). 



(98) 



(99) 



We first establish the following lemma. 
Lemma — Let nf'^ and nf'^ be the eigenvalues of p^'^^ and 



r]^'^\ respectively, and let us assume that n^^ > > 



(2) 



> • • • and n};^ > > fi)^^ > 



^2 

i = 0, 1, . . . , the following inequality holds: 

< 1. 



(2) -(2) 

n) — n) 



. Then, for all 
(100) 



Proof. We notice that by using Eqs. and ([QS]) . 

7 + ^(1^1) >o, (101) 

and 

/ + >o. (102) 

It follows from the relation r/^^) - 7/(1^1) = _ ^(101) 

that 7^(2) < p(2) ^ J g^j^^ ^(2) < ^(2) ^ 

Let A and B be arbitrary Hermitian matrices and 
let us assume that their eigenvalues are ordered in 
such a manner that > ai > • • • and > hi > • • • ^ 
respectively. If A > 5, then > hi for all i. Then, 
we have nf'^ < nf'^ + 1 and rSp < nf^ + 1. Thus, the 
Lemma is proved. (Q.E.D.) 



(2) 

Then, in the normal phase, nl = o(A^), whereas in 
the BEC phase 



.(2) 



: XiN^0{N), 



{0<Xi< 1). (103) 



B. Criterion of ODLRO for a two-component 
Fermi system in terms of irreducible graphs 

The criterion of ODLRO based on the Lee- Yang 
cluster expansion can be obtained for two-component 
Fermi systems, following procedures similar to Bose 
systems. 



(a) 



h ' ^2 , 02 



A 



A ^ 




cutting two 
internal lines 



Ti, cJi r2, 02 



Ti, CTi 12, ^2 



(b) 




V 

A A 

>x 

(2) 



FIG. 8: (a) One of the conditions appearing in the defini- 
tion of a reducible contracted 2-graph. By cutting two of its 
internal lines and associating (1,2) with the incoming lines 
and (1^2^) with the outgoing lines, the entire graph can be 
separated into two disconnected contracted 2-graphs; one in- 
cluding two incoming external lines (1,2) and two outgoing 
external lines, and the other including two incoming external 
lines and two outgoing external lines (1^, 2^). (b) An example 
of reducible contracted 2-graphs. If we cut one of its internal 
lines, the entire graph can be separated into two disconnected 
parts, (c) An example of reducible contracted 2-graphs. If we 
cut two of its internal lines, the entire graph can be separated 
into two disconnected parts which satisfy the condition given 
by in (a). 



Definition. — A contracted 2-graph is called reducible 
if it satisfies any one of the following conditions: (i) By 
cutting one of its internal lines^ the entire graph can be 
separated into two disconnected parts, (ii) By cutting 
two of its internal lines and associating (1,2) with the 
incoming lines and (l',2') with the outgoing lines, the 
entire graph can be separated into two disconnected 
contracted 2-graphs; one including two incoming ex- 
ternal lines (1,2) and two outgoing external lines, and 
the other including two incoming external lines and two 
outgoing external lines (1^,2^). 

Definition. — An irreducible contracted 2-graph is a 
contracted 2-graph that is not reducible. 

The above condition (ii) and its examples are shown 
in Fig. [HI An irreducible contracted 2-graph is called a 
simple 2-diagram in Ref. [6] . 

(2) 

In terms of irreducible contracted 2-graphs, ^r/p/.r^r^ 
can be expressed as a geometric series: 



^(2) ^^(101) ^^(101)j^(2)^(101) 

-f 7^(l^l)i^(2)^(l^l)i^(2)^(l^l) 



(104) 



15 



where 

k''J\ . = [all different irreducible 

contracted 2-graplis] . (105) 

Equation ()lQ4p can be rewritten as 

^(2) ^^(101)^^(101)^(2)^(2)^ (IQg) 



(a) 



Equations ([TQ4)) . (fTOSj) and (ITQGI) are illustrated in FigM 
From Eq. (|106p . we obtain 



Thus, the following theorem is established: 



(107) 



Theorem 3. — 
The following two criteria are equivalent: 

1) The largest eigenvalue of the two-particle reduced 
density matrix p^^^ is 0{N) as N ^ oo. 

2) The smallest eigenvalue of (r/^^*^^)) ^ — K^'^^ is 
0{l/N) asTV^oo. 

The former implies ODLRO at the two-particle level and 
the latter can be obtained from the sum over the con- 
tracted graphs. 



C. Alternative form of the criterion of ODLRO of 
a two-component Fermi system 

Following a procedure similar to Sec. IVllC. we can 

(2) 

rewrite ri , , in a geometric series which is different 

'r-,^r2;rir2 " 

from Eq. (HDH): 

^(2) + ^(2) ^ |^^(2)^' ^ |^^(2) j3 ^ . . . ^ (108) 



where 



(001) 



(100) 



^(1«)1) 



rir2;rir2 * ^y'^y'^-.yiy^ ^y'-^y'^;yiY2 '^y\y'^\yiY2 
r;r^;rir2' 



and 



(109) 



(110) 



(111) 



This geometric series is illustrated in Fig. [TOl 

(2) 

The matrix 77^, ^/ .^^^^ may be expanded in terms of its 

(2) 

eigenfunctions Xi(ri,r2) with eigenvalues as 

^ri!.i;r.r. = E ^ X* (ri , r2)x. (r'^ , r'^) . (112) 



A 



a + 




(b) 




y + t 




FIG. 9: (a) Expression of (1', 2'|ry(^) |1, 2) in terms of irre- 
ducible contracted 2-graphs in the form of a geometric series, 
(b) Expression of (l', 2'|K(2^ |1, 2) as the sum over all different 
irreducible contracted 2-graphs. 



[Note that Xi(ri,r2) 7^ Xi(ri,r2).] Since matrices r^'^^^ 
and iC^^^ commute, they can be diagonalized simultane- 
ously. Thus, we can write 



r'^r2;rir2 



^Af x*(ri,r2)Xi(r'i,r'2). 

i 

We may write Eq. ([76|) as 

00 , _-, 



(113) 



^=0 



Then, we obtain 



1 

7^' 



(114) 



(115) 



Thus, the criterion of ODLRO at the two-particle level is 
stated as the largest eigenvalue of K^'^^ becoming unity 
diS N ^ 00. 
Theorem 4- — 
The following two criteria are equivalent: 

1) The largest eigenvalue of the one-particle reduced 
density matrix p^^^ is 0{N) as ^ 00. 

2) The largest eigenvalue of the matrix K^^^ is 
1 + 0(1/7V) as TV ^00. 
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(a) 




FIG. 10: (a) Expression of (l', 2'|7y^^^|l, 2) as the sum over all 
different primary 2-graphs in the form of a geometric series, 
(b) Expression of {V ,2'\K^'^^\l,2) . 



The former implies ODLRO at the two-particle level and 
the latter can be obtained from the sum over the con- 
tracted 2-graphs. Thus we can use Theorems 3 and 4 
to find out whether the Fermi system exhibits ODLRO 
at the two-particle level based on the Lee- Yang cluster 
expansion method. 



D. Application to a uniform two-component Fermi 
gas in the tightly bound limit 

As an example of Theorem 4, we consider a uniform 
two-component Fermi gas in the tightly bound limit, 
which was discussed in Example 4 of Sec. V, and evaluate 
the critical value Zc for the Bose-Einstein phase transition 
of diatomic molecules. Using Eq. (|67|) . we have 

= z2(k;,k'2|Ti^'^)|ki,k2> + e^l^-l X (O(^) + 0{z)) 

= z2e-''K^/(4™)<5K,K'Vb(k')€(k) e^\^-\ 

+ ^2g/3|£;bix (o(^)+0(z)), (116) 

where a is an s-wave scattering length between up-spin 
and down-spin atoms. Here in the tightly bound limit, 
the following inequalities hold: < 2; ^ 1 and < 
a/ A <€. 1. According to Theorem 4, the critical fugacity 



Zc satisfies the following equation: 

l = zle''^^^\x (l + 0(^j)+0{z)y (117) 
Hence, we obtain 

e-'^l^-l/^x (l + o(^)). (118) 

This result agrees with that obtained in Example 4 of 
Sec. [Vlor in Ref [9|. The formulation given in this sec- 
tion includes the method in Ref [9] and enables us to 
systematically compute the criteria of ODLRO at the 
two-particle level. 



VIII. SUMMARY 

We have discussed the Lee- Yang cluster expansion 
method and its application to the criteria of ODLRO 
in Bose and Fermi systems. Lee and Yang [3] and de Do- 
minicis [6|] discussed graphical expansions of the grand 
partition function and the one- and two-particle reduced 
density matrices for the case of a one-component uniform 
system. We have generalized these results for the case 
of multi-component uniform and trapped systems. We 
have explained a physical meaning of primary 0-graphs. 
In particular, we have elucidated a physical meaning of 
an infinite series of the ladder-type Lee- Yang primary 
graphs, the sum of which leads to the BEG of the di- 
atomic molecules below a critical temperature. We have 
given the criterion of ODLRO of Bose systems at the 
one-particle level based on the Lee- Yang cluster expan- 
sion method. It is expressed as a converging infinite series 
of irreducible contracted 1-graphs or irreducible primary 
1-graphs. Applications to a Bose gas of hard spheres have 
also been made. Furthermore, we have given the crite- 
rion of ODLRO of Fermi systems at the two-particle level 
based on the Lee- Yang cluster expansion method. This 
is expressed as a converging infinite series of irreducible 
contracted 2-graphs or irreducible primary 2-graphs. Ap- 
plications to a two-component Fermi gas in the tightly 
bound limit have also been made. Finally, we note that 
the formulation given in this paper includes the method 
in Ref. (3] for Bose systems and in Ref. [9[ for Fermi sys- 
tems, enabling us to systematically compute the criteria 
of ODLRO. 
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Appendix A: Derivations of quantum cluster 
expansion formulae 

In this Appendix we derive the quantum cluster expan- 
sion formulae of the grand partition function ^ and the 
one- and two-particle reduced density matrices (p!3|) , (p!5]) 

and (p!7|) . We first introduce partial traces of W^^'^'^^ 

and Ui''''''^\ which are denoted as 

{yt|xi"-"^)|xt> 



1 1 



(n-f - 1)! n;! 
X (y, l^,...,n^ - ...,n;| 

Xi"^'"-^^ |x, l^., . . . , nt - 1; I4,, . . . , ni), 



(yi t,y2tl^i"^'"'^|xi t,X2t) 

(nt-2)!n;!i V ^ 

X (yi,y2,lt,---,nt-2;l;,...,n;| 
j^^"T.»4.) |xi, X2, It, . . . , nt - 2; I4., . . . , n;), 



(Al) 



(A2) 



(yi t,y2i|x("^'"^^|xi t,x2;) 

(nt - 1)! K - 1)! i ^ (A3) 

X (yi, 1^,.. .,nt - l;y2, 14,, . . . , - 1| 
j^^nT."4.) It, . . . , nt - 1; X2, 1;, . . . , - 1), 

where xi"-^'"^) = wi"""'^ or [/^"-^'"^^ and X^"^'"^) = 
Wi"''"'^ or C7^"1^'"^\ 

1. Grand partition function: Proof of Eq. ([5]) 

To prove Eq. ([5|), we first note that a general term in 
Eq. ([3]) for wl^^'^^^ is expressed as the sum over prod- 
ucts of mi^^i^{= 0, 1, 2, ... ) cluster functions Ua^'^^^ : 

(It, . . . ,iVt; • • . ,iV;|w^i'^^''^^^|it,- • ■,Nr, i^,. . .,n^) 



2,0) 



wi,o 

x(c7^''^)...C/i'-'^)) 



(A4) 



where the set of integers {ttii^^i^} satisfies the conditions 

(A5) 

(A6) 



1^=0 1^=0 



In Eq. ()A4[) . X^jm^ i } denotes the sum over all sets 
{tui^^i^} satisfying the conditions (|A5[) and (|A6[) . and 
Xlper is ^^^^ different ways of assigning 

{ri^,...,rAr^,ri;,...,rAr;} to each Ua^'^^\ In accor- 
dance with Eq. dll, we integrate wi^^'^^^ over all the 
coordinates. 

Integrating wl^^'^^^ using Eq. (|A4[) . we obtain the 



same result for every term in the sum Xlper fixed 
{^^t'Ui* "^^^ number of such terms in the sum Xlpg^. is 



n U'^wimmr 

_Zt=0 1^=0 



1 -1 



• (A7) 



By the definition of B^^^'^^^ in Eq. (j6j), we have 

We thus obtain 

W-i^^'^^)|lt,...,iVt;lt,...,7Vt> (A9) 

= (NMN.i) E n n ^ {Bi^^^^^f'^''^ . 

The grand partition function can be obtained by substi- 
tuting Eq. (|A9l) in Eq. (g]): 



oo oo 



= =E E E 

oo CX) CX) ^ 

= nn E ;;nT7(4'4'f """•)'""" 



CX) CX) 



=exp EE4^4^^^'^^'^M- 



(AlO) 
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2. One-particle reduced density matrix: Proofs of 
Eqs. ([T3|) and ([T4]) 

To prove Eq. ([13]), we first rewrite tiie one-particle re- 
duced density matrix using Wi^^'^^\ We note that 



E 



(N^) (JVi) 



X (tA.|*|(x)|l,...,Ar-l) 

X {l,...,iV-l|*t(y)l^.)e"''^'*' 



= iV,. d^^-^-^^r, d^-^r, (All) 

x^/>*(x, lt,...,iVt-l;li,...,iVi) 



(N) 



X (y,l^,...,iV^-l;l;,...,iVJ 

wi'''''''^\x,M,...,N^-l;h,...,N^). 

Then, the one-particle reduced density matrix is rewrit- 
ten as 

(xt|p(^)|yt) := {*|(x)*t(y)) 

Using the definitions of the cluster functions (|3]), 

Wi^'^^'^ in Eq. jSn]) can be expressed as the sum over 
products of the cluster functions: 

(y,l^,...,7V^-l;l;,...,7V;| 

iyi^^'^^^|x,l^,...,7V^-l;l;,...,7V;) 

= E E E E 

n^ = l n^=0 {mz^,z^} per. 

(y, ai, . . . , ttn^-i; ai, . . . , an^| 

y-y-(n-|-,n^) I \ 

Ua |x, ai, . . . , ttn^-i; ai, . . . , an^/ 
x(6i|C/(i>°)|6i)---(5™,JC/(i-°)|5^,.„) 

X • • • 

X (ci, . . . , Q^; ci, . . . , cijui^^'^^^\ci, . . . , Q^; ci, . . . , q^) 



(A13) 



where ^ 



is the sum over all sets of integers 



{mi^^ij satisfying = + ^^^h,h^ and ^ 



is the sum over different ways of grouping as follows: 
(ai, . . . , an^_i), . . . , (6^i,o)}^ • • • 



{(ci,...,Qj,...,(...,Q^m,^,,J},...: 



(A14) 



where ai, . . . , an^-i, ... is a permutation of the coor- 
dinates {1^, . . . , A^-^ — 1}, and 



(ai, . . . , an^-i), . . . 

{(Ci, . . . , C^^), ...,(..., U-^mi^j^)}j • • • •) 



(A15) 



where ai, . . . , an-^^-i, ... is a permutation of the coordi- 
nates {1;, . . . , A^;}. In the groupings (|A14p and 
within each pair of round brackets the integers are ar- 
ranged in ascending order and within each curly bracket 
the round brackets are arranged such that their first en- 
tries follow ascending order. In general, Eq. (|A13[) is 

a product of one U^^'^^^ which includes x and y as 
two of its variables, and products of cluster functions 
Ua^'^^\ which include other particles as variables. Ac- 



cording to Eq. (|A12[) . we need to integrate Wa 
Using Eq. (|A13p . if we integrate Wa 



over coordi- 



nates {1^, . . . , — 1; 1|, . . . , A'l}, we obtain the same 
result for every term in the sum ^^^^ for fixed n^, 
and {mi^^i^}. The number of such terms in the sum Xlpg^ 
is 



(n^ - l)!n;! 



We thus obtain 

(ytlW^r^-'^'^Ixt) 

= E E E (ytP^ 

- n n ^ (^^'^''^0 

Substituting Eq. (|XT6]) into Eq. ([XT2]) . we obtain 

(xtlP^^^lyt) 



(A16) 



lE E E E E 



|xt) 



(A17) 



Using Eq. (|XTO|) . we find 
(xtlp('Vt) 



CO oo 



per. 



= E E4'"\"'(ytit^"^'"^^ixt). 

The proof of Eq. ([Tlj) can similarly be made. 



(A18) 
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3. Two-particle reduced density matrix: Proofs of 

Eqs. ([iiD-idzD 

To prove Eqs. (p!5 ]) - (pT|) . we rewrite the two-particle re- 
duced density matrix using in a manner similar 
to what we have done to the one-particle reduced density 
matrix in the preceding subsection. The results are 



(xi t,x2t Ip^^^Iyi t,y2 1) 
:= (^|(xi)^|(x2)^t(y2)^t(yi)) 



and 



^ OO CX) 

X (yi t,y2t|W-i'^'''^'^|xi t,x2t>, 



(Xi t,X2; 1/9(2) |yi t,y2 i) 

:= (M(xi)*Ux2)*;(y2)*t(yi)) 



(A19) 



-J CXD OO 

E E E 4' 



A/'t-l AT^-l 



x(yi t,y2;|V^I^^'^^^|xi t,x2;). 



(A20) 



First, we prove Eq. (p!5]) . We now use the definitions 
of the cluster functions ([3]) and note that each wi^^'^'^^ 
in Eq. (|A19[) can be expressed as the sum over products 
of the two forms as follows: 

{yi,y2,ai, ... ,an^-2;bi, ... ,bnj 

^K'^i) |xi, X2, ai, . . . , an^_2; 6i, . . . , b^^) 



In Eq. (|A21[) . the set of integers {tui^^i^} satisfies A^cr = 
Tier + ^T'^i^cr^h^h^ (E22]) the set of in- 

tegers {mi^^i^} satisfies the conditions Ncr = nia + 
^2(7 + YlT^^i^^'^h.h' general, (i) the term (|A2ip 
is a product of one U^^'^^^ which includes xi, X2, yi 
and y2 as four of its variables, and products of clus- 
ter functions Ua^'^^\ which include other particles as 
variables, and (ii) the term ()A22[) is a product of one 
^^ni^,ni^) ^]^^^^]^ includes xi and yi as two of its vari- 
ables, one /7^^2t»^2i) includes X2 and y2 as two of 
its variables and products of cluster functions Ua^'^^\ 
which include other particles as variables. In accor- 
dance with Eq. (|A19|) . we integrate W^^'^^\ Using 

Eqs. (|A2T]) and (|A22]) . we integrate wi^^'^^^ over co- 
ordinates {1^, . . . , — 2; 1|, . . . , N^}. We thus obtain 

(yi t,y2t|VFi'^^''^^^|xi t,x2t) 

= T.T. E<ylt,y2t|c/^'"^^|xlt,x2t) 
+ E E E E E 

(yi t |l7^^'"^^^|xi t)(y2 t |i7("^^'"^^^|x2 t) 



(A23) 



X (di, . . . , ; ei, . . . , ei \Ua^'^^^ . . . , ; ei, . . . , ei ) 



•••(..., di^mi^j^ 5 • • • 7 ^limi^j^ I 



I • • • 7 di^rni^j^ 5 • • • 7 ^limi^j^ ) 



X • 

and 



(A21) 



(yi 

|xi, ai, . . . , ttni^- 
X (y2, Ci, . . . , Cn2-|~-i; • . • , <^n2^ I 

x(ei|C7(i>°)|ei)---(e™,.JC7(i>°)|e^,.„) 



(A22) 



Substituting Eqs. (lAlOl) and (lA23]) into Eq. (|AT9]) . we 

have 

(xi t,x2t|p('^|yi t,y2 1> 

OO OO 

= E E 4~'h'(y^ t,y2 1 |c7^'"^)|xl t,x2 1> 

/ CO CO \ 

+ E E 4^'"'C'(yi^i^^"""^i^it> 

/ OO OO \ 

X E E <""\""(y2tiol"^^'"^^)|x2t> . 

(A24) 



Then, using Eq. (jAlSp we obtain 

(xi t,x2t|p('^|yi t,y2 1) 

= (xitb'')|yi t>(x2t|p(')|y2 1) 



(A25) 



OO OO 



J2 E 4~'h'(y^ t,y2 1 iol"-"^)|xi t,x2 1>. 
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The proof of Eq. (p!6]) can simiarly be made. 

Next, we prove Eq. (fT7|) . Following a procedure similar 
to the above derivation of Eq. ([15]) , we obtain 

(yi t,y2i|wi'^^''^^^|xi t,x2i) 

= E E E t,y2 i |/7^'"^)|xi t,x2 ;) 

A/'t AT, 

^ n n 

+ E E E E E 

(yi t |c7^^^'"^^^|xi t){y2 i \ui"'""''^\^2i) 

A/'t AT^ 

X n T]^—(b(^^mY'''\ (A26) 



Substituting Eqs. (|XTO|) and (|X26ll into Eq. (iX20|) . we 
find 



(xi t,x2;|p(')|yi t,y2 ;) 



CX) oo 

= E E <'"'<'"'(yi t,y2 i ic/r^-^^|xi t,x2 ;) 

/ CX) CX) 



E E <^^"'C^(yit|c/r^'"^^^|xit) 

, ni^ = l ni^=0 



CX) CX) 



E E <''<''"'(y2i|f^i"'''"''^|x2i) . 



V n2t=0 n2 I =1 



(A27) 



Then, using Eq. ()A18|) we have 

= (xit|p('^|yi t)(x2;|p«|y2 ;) 

CX) CX) 



(A28) 



^ ^ 4^-^.;^-^yi t,y2 i t,x2 



4. Mathematical justification of quantum cluster 
expansion 

The above derivations can actually be mathematically 
justified. We consider a uniform system or a trapped 
system. To be rigorous, let us assume that the system 
is confined in an infinite potential well with a volume V 
and the interaction v has a hard repulsive core. Then 
^{N^^N^) yg^j^^gj^gg sufficiently large values of or 
A^^ and we denote their upper bounds as and M^, 
respectively. 



First, we demonstrate the validity of Eq. (j5j). The 
grand partition function reads 



Ml 

^ Nfl Nil 



)(A^t,A^i) 



(A29) 



where 

Q(^t.A^.) I d3iVtr, I d'^^r.m, . . . , iV,; 1,, . . . , N,\ 

wi'''''''^m,...,Nr,U,...,N^), (A30) 

is a partition function. Equation (jA29p is positive and is 
thus a polynomial of and with no zeros on both the 
positive real axes. Therefore, logS(2;^, z^) is holomorphic 
near the origin and along the positive real axes in both 
the complex z^ plane and z^ plane. Near the origin this 
logarithm can be expanded as a power series, and this 
power series is nothing but Eq. ([5]). If we understand 
Eq. ([5]) to mean the analytic continuation, Eq. (j5j) is valid 
for all positive values of z^ and z^. 

Next, we show the validity of Eq. Since W^^^'^^^ 
vanishes for sufficiently large values of or and S is 
a polynomial function of z^ and z^ with no zeros on their 
positive real axes, from Eq. ()A12p the one-particle re- 
duced density matrix (x t |y t) is a rational function 
of z^ and z^. Therefore, the one-particle reduced density 
matrix is holomorphic near the origin and along the pos- 
itive real axes in both the complex z^ and z^ planes. 
Near the origin the one-particle reduced density matrix 
can be expanded as a power series, and this power se- 
ries is exactly Eq. ([T3|) . If we understand Eq. ([13]) to 
mean the analytic continuation, Eq. (p!3|) is valid for all 
positive values of z^ and z^. The two- and multi-particle 
reduced density matrices can be mathematically justiffed 
in a similar manner. 

Note that, for the above justification, interaction v 
needs to have a hard-repulsive core, but the 5- wave scat- 
tering length can be either positive or negative. 

Appendix B: wL^^'^^^ in terms of Vl^^^t'^^^ 

We prove Eq. ([2T]) . The eigenfunctions of 
}{iN^,Ni) classified according to the ten- 

sor product of the irreducible representations of 
the permutation group of objects and that 

of TV; objects. If . . . , A/"!; 1^, . . . , 7V|) be- 

longs to an irreducible representation (g) D^, 

then e^e^V^,(P(l^),...,P(A^|);Q(l^),...,Q(7V|)) also 
belongs to the same representation (g) D^. Hence 

^ J2 e^e«^i(P(4), . . . , P{N{y, Q{1[), Q(iV|)) 

PeSN^QeSN^ 

(Bl) 

belongs to ^ D^. However, the function (jBip is sym- 
metric. Hence, if is not the symmetric repre- 
sentation, (|Bip is zero. On the other hand, if (g) is 
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the symmetric representation, then 

P€Sn^ QeSN^ 

= iN^\)iN^\)Ul\, ...,N{;1[,..., N[). (B2) 
Using the definitions ([2]) and (fTOj) . we obtain Eq. (|2T]) . 

Appendix C: Correspondence between Lee- Yang and 
our notations 

In this Appendix we hst the correspondence between 
the Lee- Yang's notation [3] and ours in momentum space. 

For simphcity, we consider a one-component system. 
The r]Q^\k.) and 7^^"^^(k) functions used in this paper are 
defined by //o^^(k) := (1 — eze~^^'^)~^ and r]^^\]i.) := 



(^(k)^''"(k)). The corresponding m(k) and M(k) func- 
tions used in the original paper of Lee and Yang (sf 
are defined by m(k) := z{l — eze~^^ and M(k) := 
z{l - (^t(k)^(k))}. Thus we obtain m(k) = zr]'^\]<i) 
and M(k) = zrf^^^i)^). These relationships are listed in 
Table III 



TABLE L Correspondence between Lee- Yang and our nota- 
tions 





Lee- Yang this paper 


solid line 
dotted line 
thick solid line 


z 1 
m(k) ry^^)(k) 
M(k) 77(^)(k) 


/-vertex 


(...|T^^^|...) (...|T^^^|---) 
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